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We report tests of local position invariance based on measurements of the ratio of the ground 
state hyperhne frequencies of 133 Cs and 87 Rb in laser-cooled atomic fountain clocks. Measure- 
ments extending over 14 years set a stringent limit to a possible variation with time of this ratio: 
dln(^R,b / Vcs) / 'dt — (—1.39 ± 0.91) x 10 -16 yr _1 . This improves by a factor of 7.7 over our previous 
report [H. Marion et al, Phys. Rev. Lett. 90, 150801 (2003)]. Our measurements also set the first 
limit to a fractional variation of the Rb/Cs frequency ratio with gravitational potential at the level 
of (?d\a{vBh/v Cs )/dU = (0.11 ± 1.04) x 10" 6 , providing a new stringent differential redshift test. 



The above limits equivalently apply to the fractional variation of the quantity a 



x (flRb/ffCs 



which involves the fine-structure constant a and the ratio of the nuclear g-factors of the two alkalis. 
The link with variations of the light quark mass is also presented together with a global analysis 
combining with other available highly accurate clock comparisons. 

PACS numbers: 06. 30. Ft, 06.20.Jr, 04.80.Cc 



Einstein's equivalence principle is one of the founding 
principle of general relativity. Many experiments have 
been dedicated to testing the validity of this principle [l[ , 
several of them searching for variations of fundamental 
constants, cither on cosmological timcscalcs using astro- 
nomical and gcochcmical data or in our present epoch, 
by exploiting highly accurate atomic clocks. Such varia- 
tions would violate local position invariance (LPI) , one of 
the three components of Einstein's equivalence principle. 
The possibility that dimensionless fundamental constants 
might change in time or space is allowed or predicted by 
alternative theories aimed at unifying gravitation with 
the other fundamental interactions, hence the strong in- 
terest in this search, which could reveal physics beyond 
general relativity and the Standard Model of particle 
physics. There has been a claim of a variation of the fine- 
structure constant a over cosmological time scales from 
the analysis of Quasar absorption spectra 2j. This was 
followed by other conflicting observations 13 1 and by sev- 
eral controversies (see, for instance, Ref. [J]). Recently, 
it was suggested that conflicting observations could be 
reconciled if one assumes that a varies in space rather 
than in time [||. Atomic clocks offer the possibility to 
search for variation of constants at the present epoch in 
laboratory-based experiments whose interpretations are 
fully independent of any cosmological model 

In this letter, we present new LPI tests obtained 
by comparing the ground state hyperhne frequencies of 
133 Cs and 87 Rb atoms over a period of ~ 14 years. Our 
measurements give the second most stringent limit to 
date to a possible variation with time of the ratio of two 
atomic frequencies, and the most stringent for two hy- 
pcrfinc frequencies. They also provide the first limit to a 
variation of the hyperhne frequency ratio ^Rb/^Cs with 
gravitational potential. The link with fundamental con- 



stants of the Standard Model is made with atomic and 
nuclear structure calculations (l6l . 17 1. Combining our 



measurements with other clock experiments, we set up- 
dated limits to variations of the fine-structure constant 
a, of the electron-to-proton mass ratio fj, = m e /m p and 
of the ratio of the light quark mass to the Quantum Chro- 
modynamics mass scale m q / Aqcd- 

Our measurements exploit the LNE-SYRTE atomic 
fountain ensemble schematized in Fig. [TJ which consists 
of three 133 Cs Primary Frequency Standards FOl, F02- 
Cs and FOM, sharing a common ultra-low noise cryo- 
genic oscillator. Among them, F02 is a dual fountain 
that can also operate with 87 Rb. The three Cs clocks 
operate on the hyperhne clock transition at 9.2 GHz and 
the Rb clock on the corresponding 6.8 GHz transition. 
A detailed description of the latest developments of this 
fountain ensemble is given in Ref. [HI . Notable advances 
relevant to the present work were the simultaneous op- 
eration with Rb and Cs in F02 [l!| and large improve- 
ments in reliable, unattended operation for all fountains, 
allowing for quasicontinuous Rb/Cs comparisons for sev- 
eral months. Also, fountain accuracy improved along the 
years, notably with the several recent studies of system- 
atic shifts [IBS]. In F02, Rb and Cs atoms are si- 
multaneously laser-cooled, launched, state-selected, and 
probed with the Ramsey interrogation method, and fi- 
nally selectively detected by time-resolved laser-induced 
fluorescence, in the same vacuum chamber [l||, as shown 
in Fig. [1] Rb and Cs are launched at slightly different 
velocities to separate the two clouds during interroga- 
tion and thereby avoid interspecies collisions. Typical 
accuracies over the recent years relevant to this work 
are 4 - 5 x 1CT 16 for FOl, F02-Cs and F02-Rb, and 
7-8 x 1CT 16 for the transportable fountain FOM. The de- 
tailed systematic uncertainty budgets for the four clocks 
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FIG. 1. LNE-SYRTE atomic clock ensemble. The block di- 
agram shows the 3 fountains FOl, F02 and FOM connected 
to the common ultra-low noise cryogenic sapphire oscillator. 
The F02 fountain, schematized on the right, can operate with 
S7 Rb and 133 Cs simultaneously, enabling the present hyperfine 
frequency comparisons. The optical clocks, the H-masers, the 
sapphire oscillator, as well as the links to other remote labora- 
tories and to international timescales (TAI is Temps Atomique 
International) via satellite time transfers (GPS, TWSTFT) 
provide for further tests of Einstein's equivalence principle 
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are available in Table III of Ref. [3]. The four clocks 
measure the frequency of the common local oscillator 
(Fig. [I}, with cycle times ranging from 1.1 s to 1.6 s. 
The data from each clock are corrected for all systematic 
shifts (as described in Ref. [l8j], Sect. IV. A) and aver- 
aged over synchronous intervals of 864 s. Next, for each 
available Cs fountain, the Rb/Cs frequency ratio is de- 
termined over these synchronous intervals, which removes 
the frequency of the common local oscillator. Typically, 
the fractional frequency instability of a Rb/Cs compar- 
ison is limited by quantum projection noise, and ranges 
from 8 x 10 -14 to 1.2 x 10~ 13 at 1 s. A fractional reso- 
lution of 2 x 10 -16 for the t'Rb/^'Cs ratio is reached in a 
few days of averaging time, after which the overall uncer- 
tainty of the comparison becomes limited by systematic 
uncertainties. 

Figure [2] presents the temporal record of the Rb/Cs 
comparisons since the first operation of F02-Rb in 1998. 
Up to 2008 (the first six points), F02 was operated ei- 
ther with 87 Rb only or alternately with 87 Rb and 133 Cs 
atoms, and the Cs reference was FOl or FOM, or both. 
Since the end of 2008, F02 has been operating with Rb 
and Cs simultaneously and the dominant Cs reference is 
F02-Cs. Each point in Fig. [5] represents the average over 
a duration ranging from a few weeks to several months, 
the averaging periods being chosen based on significant 
modifications on either fountain. When two or three Cs 



FIG. 2. Temporal record of fractional variations of the 
i / Rb/i / Cs hyperfine frequency ratio. The error bars are the 
total 1 a uncertainties, dominated by the systematic uncer- 
tainties. The horizontal bars show the duration of each com- 
parison. The solid line is the weighted fit to a line with 
inverse quadratic weighting. The origin of the vertical axis 
corresponds to the S7 Rb secondary representation of the SI 
second recommended by BIPM, with a recommended uncer- 
tainty 3 x 10" 15 (dashed lines) 

references are available, we compute a weighted average 
taking into account the total uncertainty and the amount 
of measurement time for each Cs/Rb pair. Error bars in 
Fig.[2]are the overall one standard deviation (1 a) uncer- 
tainties, which are dominated by systematic uncertain- 
ties. 

The solid red line in Fig. [5] is the result of the weighted 
linear least-squares fit of a straight line to the data with 
inverse quadratic weighting, i.e. weights inversely pro- 
portional to the square of the error bars: 

1 in ( ^ = (-1.36 ± 0.91) x 10- le yr -1 - (1) 
at V ^Cs / 

The yearly drift deviates from zero by 1.5 er which is not 
statistically significant. The uncertainty of this result 
improves by a factor 7.7 over our 2003 report 0- It is 
one of the most constraining accurate clock comparisons 
over long timescales (see Table ITT]) . 

Hyperfine splitting (hfs) energies scaled to the Rydberg 
energy (Roohc) depend on the fine-structure constant a 
and on the nuclear g-factor. Consequently, the VRb/vcs 
hyperfine frequency ratio is sensitive to variations of 
a and of the nuclear g- factors of the two atoms [(§. 
From the recent atomic structure calculations of the a- 
dependent relativistic effects [l(|, we get dln(i/Rj :) / 1/ Cs) = 
dln[ a ~ 0A9 x (W5Cs)]- Thus, Eq.fflyields 

-In ( a- 49 x ^) = (-1.36±0.91)xl0- 18 yr" 1 - (2) 
at V SCs / 
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FIG. 3. Search for a modulation of the vwd/vcs hyperfine fre- 
quency ratio synchronous with the annual change of the solar 
gravitational potential on Earth. The red curve is the fitted 
modulation magnified by 10 for visibility. The inset shows the 
fitted amplitudes in phase (X) and in quadrature (Y) with 
the potential modulation. Fits with fixed phase (red dots) 
or arbitrary phase (green triangles) , with (filled symbols) and 
without (unfilled symbols) accounting for a slope. 



The link between the (/-factors and the fundamental pa 



rameters of the Standard Model was achieved in Ref. [17j , 
where it is shown that a variation in g-factors can be re- 
lated to a variation of the light quark mass m q scaled to 
the Quantum Chromodynamics mass scale Aqcd- The 
relation for the Rb and Cs ^-factors yields 



|ln(a- 049 (/^; 
at Aqcd 



-0.021) = (_i. 36± o.9i)xlO- 16 --- 1 



(3) 

showing that hyperfine transitions can test both the elec- 
troweak (a) and the strong (?ti 9 /Aqcd) interactions, al- 
though without distinguishing between the two contribu- 
tions. 

Next, we search for a possible coupling of the VRb/^Cs 
frequency ratio to gravity. For that purpose, we group 
our Rb/Cs data over shorter intervals of typically one 
month, as shown in Fig. [31 and we look for variations 
correlated with the annual change of the Sun gravita- 
tional potential on Earth, AU(t) « AU cos[f2 ffi (i — t p )], 
where Q,® is the Earth's orbital angular frequency, t p 
a perihelion date, and AU s» GM Q e/a, to first or- 
der in the Earth's orbit eccentricity e = 0.0167. G = 
6.674 x 10 _11 m 3 kg -1 s -2 is the gravitational constant; 
Mq = 1.98 x 10 30 kg the solar mass; and a = 1 AU the 
orbit semi-axis, yielding AU/c 2 ~ 1.65 x 10~ 10 , where 
c is the speed of light. We use the sign convention 
recommended by the International Astronomical Union 
prjj ; i.e., U(r) +GM@/r is positive and maximum 
at perihelion. The fit function to the data in Fig. [3] is 



cos [1^0 (t — t p )], the modulation amplitude being a free 
parameter. For t p , we take the Jan 3 rd , 2010, perihelion: 
MJD=55199. The fitted amplitude is (0.18±1.71) x 10~ 16 
for the fractional modulation of ^Rb/^Cs ratio, which 
yields 



c 2 d ln (^_\ (Q n ± LQ4) x 1Q -6 



(4) 



This result can be interpreted as a differential redshift 
experiment, testing whether the gravitational redshift 
depends on the clock composition. The degree of LPI 
violation is usually parameterized by a composition- 
dependent coefficient j3 in the gravitational redshift [l[ 
dv/v = (1 + P)dU/c 2 . Our measurement yields 



/3( 87 Rb) - p( li6 Cs) = (0.11 ± 1.04) x 10" 



(5) 



The previous most sensitive test of this type is based 
on comparisons between H-masers and Cs fountains [l2| ■ 
Our result is 1.4 times more stringent. Further, it in- 
troduces a new atomic transition ( 87 Rb hfs). Finally, it 
relies on clocks with high accuracy for both Rb and Cs, 
while H-masers are known to be influenced by "poorly 
understood and unmodeled factors" fl2| . 

A violation of LPI in this test could occur via a 
coupling of fundamental constants to gravity, as pre- 
dicted by some unification theories. Using the pre- 
viously mentioned sensitivity of the vnb/^Cs ratio to 
variations of a and of the nuclear ^-factors, we obtain 
c 2 dln(a- a49 x.g Rb / 5Cs )/d(7= (0.11 ± 1.04) x lO" 6 , where 
again the variation of <7Rb/<?Cs can also be presented as a 
variation of (m q /AQCD) -0 ' 021 - 

Our experiment can be exploited to bound a possi- 
ble variation of constants in space. In the Barycentric 
Celestial Reference Frame, the Earth approximately fol- 
lows a circular orbit of radius a — 1 AU at the angular 
frequency £1$. To search for a variation of Ush/^Cs corre- 
lated with this motion, we now fit the data of Fig. [3Jwith 
cos(f2®i + tp) with an arbitrary phase tp. We find a mod- 
ulation amplitude A ln(^ Rb /^c s ) = (0.23 ±1.80) x 10~ 16 . 
The inset of Fig. [3J shows the result of this fit in the form 
of in-phase (X) and in-quadrature (Y) components with 
respect to the perihelion. From this analysis, we find 
that dlnfab/ vcs) /dr = dln[a~ - 49 x (g Rh /g Cs )]/dr = 
(0.23±1.80) x 10" 16 AU- 1 or (0.15 ±1.20) x 10~ 27 m" 1 . 
Another, probably more relevant, approach is to consider 
the motion of our experiment with respect of the rest 
frame of the cosmic microwave background, as suggested 



in 27j. This motion is largely dominated by the motion 
of the solar system at a velocity of 369 ± 0.9 km.s -1 
in the direction (168°, -7°) [28[. In the cosmic mi- 
crowave background rest frame, the result of Eq.Q] leads 
to dln(v Rh /v Cs )/dR = dln(a" ' 49 x g Rh /g Cs )/dR = 
( — 1.17±0.78) x 10~ 29 m _1 , which is still approximately 
two orders of magnitude short of resolving the spatial 
variation suggested in Refs. [El l27j]. 
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TABLE I. Sensitivity coefficients k a , k U; k ? of atomic transition frequencies used in current atomic clocks to a variation of a 
[23l. [2^]. of — m e /m p and of ?ti 9 /Aqcd [I6lll7l|. These transitions are hyperfine transitions for 1 Hhf s , S7 Rb, 133 Cs, and optical 
transitions for 1 H(1S — 2S) and all others excepted Dy. For Dy, the rf transition between two closely degenerated electronic 
levels of opposite parity is used in the two 162 and 163 isotopes [ToL [Til . [25j . 





87 Rb 


133 Cs 


Hhfs 


^(lS^S) 


171 Yb+ 


199 Hg + 


87 Sr ( 162 Dy- 163 Dy) 


27 A1+ 


hot 


2.34 


2.83 


2.0 


~0 


1.0 


-2.94 


0.06 1.72xl0 7 


0.008 




1 


1 


1 

















kg 


-0.019 


0.002 


-0.100 


















TABLE II. Results of atomic frequency ratio measurements used to search for variations of fundamental constants. Second 
column: sensitivity coefficients deduced from Table U for each atomic frequency ratio X. Third column: variation with time in 
yr _1 . Last column: variation with the gravitational potential. We give the measured fractional modulation of the frequency 
ratio X scaled, to solar gravitational potential modulation AU/c 2 ~ GMq£/ '(ac 2 ) « 1.65 x 10 -10 . 



Frequency ratio X 






hq 


d\n{X)/dt (yr- 1 ) 


Ref. 


c 2 d\n(X)/dU 


Ref. 


Rb/Cs 


-0.49 





-0.021 


(-1.36 ±0.91) x 10" 16 


This work 


(0.11 ± 1.04) x 10" 6 


This work 


H hfs /Cs 


-0.83 





-0.102 






|0.1 ± 1.40 x 10~ G 


[12] 


H(1S - 2S)/Cs 


-2.83 


-1 


-0.002 


(-32 ±63) x 10~ 16 


M 






Yb+/Cs 


-1.83 


-1 


-0.002 


(-4.9 ±4.1) x 10~ 16 


[15] 






Hg+/Cs 


-5.77 


-1 


-0.002 


(3.7 ±3.9) x 10~ 16 


m 


(2.0 ±3.5) x 10^ 6 


® 


Sr/Cs 


-2.77 


-1 


-0.002 


(-10 ±18) x 10~ 16 


[13J 


(-11.5 ± 18.2) x 10 -6 


[13] 


( 162 Dy- 163 Dy)/Cs 


1.72xl0 7 


-1 


-0.002 


(-4.0 ±4.1) x 10 _s 




(134 ± 104) 




Al+/Hg+ 


2.95 








(-0.53 ±0.79) x 10" 16 


M 







A possible fractional variation of any atomic frequency 
X can be related to variations of three dimensionlcss 
constants, a, li = m e /m p , and m q /AQcr>. We can 
write <im(X) = fc Q dln(a)±fc M dln(/i)±fc 9 (iln(mg/AQCD), 
where the fc's represent the sensitivity coefficients of the 
specific transition X. This offers the possibility to set in- 
dependent limits for the three constants using multiple 
atomic transitions. Next, we provide such a global analy- 
sis combining our measurements with other reported ac- 
curate clock comparisons. Table U lists the sensitivity co- 
efficients that we use in this analysis and which are taken 
from the most recent atomic and nuclear structure cal- 
culations 1|1 17, 23, 24 1 . Sensitivity to li and m g /AQCD 
comes from the nuclear magnetic moment involved in 
hyperfine transitions. Both optical and hyperfine transi- 
tions have dependence in the fine structure constant a via 
relativistic corrections. Table ITll gives the corresponding 
sensitivities for clock comparisons of relevance here. De- 
pendence in li cancels out in hyperfine-to-hyperfine com- 
parisons. Optical-to-optical comparisons are sensitive to 
a alone, as exemplified here by the Al + /Hg + comparison. 

To set independent limits to time variations of con- 
stants, we perform a weighted least-squares fit to all ex- 
perimental results listed in Table ITT1 (3 rd column), includ- 
ing our result of Eq. [TJ This fit yields independent con- 
straints for the three constants a, li and m g /AQCD as 
reported in the first row of Table IIIII The constraint rel- 
ative to a is mainly determined by the Al + /Hg + compar- 
ison. In this fit, only the Rb/Cs comparison disentangles 
(i and 77j,j/Aqcd- It is therefore essential to constrain 



m q / Aqcd- This stems from the fact that optical fre- 
quency measurements are all performed against Primary 
Frequency Standards, i.e. against the Cs hyperfine fre- 
quency. We note that the constraint for a is slightly less 
stringent than in Ref. 14[ because we are using the more 
recent and reduced sensitivity coefficient of Ref. [2^ . 

Similarly, we perform a global analysis for the varia- 
tion with gravitational potential exploiting the available 
comparisons (Table HH last column). We could find that 
for the Sr/Cs comparison, the modulation amplitude of 
the gravitational potential was overestimated by a fac- 
tor of 2 relative to the frequency modulation in Ref. [l3[ . 
The table gives the corrected, 2 times less stringent con- 
straint for this comparison. Also, we have checked the 
consistency between the conventions that we have cho- 
sen (sign of U, phase of the modulation) and those in 
Refs. [3, 1 1 lN 1 3 1 . The least-squares fit to these results 



yields independent constraints for the three couplings to 
gravity, reported in the second row of Table IIIII 

Finally, we note that an alternative approach consists 
in using a, m e /AQCD, and m g /AQCD instead of a, li, 
and m g /AqcD as parameters of the Standard Model. In 
this other approach, the sensitivity coefficients of Table 
U become k' a = k a , k' e = and k' q = k q - 0.048 fll| . 
The final results of this second analysis is identical for 
a and to 9 /Aqcd, and we find that dln[m e /AQCD]/^ = 
(4.9±3.7)xl0- 16 yr _1 and c 2 dln[TO e /A QCD ]/dC/ = (-4± 
17) x 10~ 6 . 

We have reported highly sensitive tests of LPI using 
87 Rb and 133 Cs atomic fountain clocks. Exploiting also 
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TABLE III. Results of the global analysis of the atomic clock 
comparisons given in Table [TTJ constraints on temporal vari- 
ations and couplings to gravitational potential for the three 
fundamental constants: a, /i = m e /m p , and m 9 /AQCD- 





ln(a) 


ln(/i) 


ln(m 9 /A QC D) 


d/dt (xl0- 16 yr _1 ) 


-0.25 ±0.26 


1.5 ±3.0 


71 ±44 


C 2 d/du (xicr 6 ) 


0.5 ±2.7 


-4± 15 


-5 ±28 



other available clock comparisons, we set stringent con- 
straints to possible variations of a, fx = m e /m p and 
m g /AQCD with time and gravitational potential. The 
rapid developments seen in atomic clocks, especially in 
the optical domain, as well as the dramatic improvement 
in long distance clock comparisons allowed by coherent 
optical fiber links 0, [3(| or the ACES space mission (3l| 
will largely diversify and enhance these tests. 
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